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Abstract 

In a previous paper pQ, we studied the rj mass and formulated its chirally symmetric coupling 
to fermions which induces electric dipole moment (EDM). Here we calculate the EDM to one- loop. 
It is finite, having no ultraviolet divergence while its infrared divergence is canceled by soft photon 
emission processes exactly as for 9 = 0. The coupling does not lead to new divergences (not present 
for sin# = 0) in soft photon processes either. Furthermore, as it was argued previously pQ, the 
EDM vanishes if suitable mixed quantum states are used. This means that in a quantum theory 
based on such mixed states, a strong bound on EDM will not necessarily lead to a strong bound 
such as I sin 9\ < 10 -11 . This fact eliminates the need to fine-tune 6 or for the axion field. 
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I. INTRODUCTION 

A non-zero electric dipole moment (EDM) of a nucleon implies parity (P-) and time- 
reversal (T-) violations. In the conventional approach to QCD, its #-term in the action, 
which violates P and T, induces an electric dipole moment d^ of the neutron. For small 8, 
it is [2H6] 

\d N \ « 9 ■ 2 • 1(T 16 cm. (1) 
On the other hand, the current experimental bound on is 

|djv| < 6 • 1(T 26 cm (90% confidence level), (2) 
suggesting that 6 is nearly zero: 

|0| < 10" 10 radians. (3) 

The need for such fine-tuning of 9 is the strong CP-problem. 

A point of this work is that d^ vanishes for any sin# if appropriate mixed states are 
used. The general considerations of our previous work [7] is thus confirmed by an explicit 
calculation. 
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As a first step to show this result, we also implement a chiral model calculation of EDM 
for the conventional QCD ^-vacuum. It is finite, having no ultraviolet divergence while its 
infrared divergence is canceled by soft photon emission processes exactly as for 9 = 0. The 
coupling does not lead to new divergences (not present for sin 9 = 0) in soft photon processes 
either. We believe that this model which fully incorporates chiral symmetry merits attention 
and can be a useful tool to investigate P- and T-violations. 

Our general one-loop result for the P- and T-violating fermion-photon vertex is reported 
in later sections. For photons on mass shell, and for the conventional pure QCD states, it 



gives the EDM (see equation (63) ahead) 

\d N \ = 8vr« M2 (4) 

for a fermion of mass M, with a = e 2 /AiT = 1/137. It depends also on a new mass scale 

In section II, we briefly recall our results from the earlier paper p.]. In particular, we 
explicitly show that the model contains the analogues of "covariant" and "consistent" axial 
charges [SI [9]. As in QCD, while the former is not conserved, the latter is. 

The Dirac equation and the fermion propagator are modified by the #-angle. They are also 
described in section II. We also need an expression showing the dependence of Dirac spinors 
on sin#. We find that as well in section II. Finally in this section, we sketch perturbation 
theory and show that the #-angle modifies either the Dirac propagators or the external 
spinors, but not both. 

Section III takes up the analysis of the one-loop diagram for EDM. It is ultraviolet finite. 
It contains an infrared divergent contribution which is canceled by soft photon emission 
processes exactly as in standard QED. The remaining 0{ot) contribution to EDM is evaluated 
to all orders in cos 9. 

In section IV, we recall the use of appropriate mixed states and how it alleviates P- and 
T-violations due to non- vanishing sin 9. It is then applied to the QCD 9- angle to show that 
EDM vanishes for these mixed states. 

In the final section V, we make general remarks in the use of mixed states for the restora- 
tion of anomalous symmetries. 
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II. RESUME AND REMARKS ON PREVIOUS WORK 



In this section, we summarize the pertinent results from our previous work [T] on the r\' 
effective Lagrangian. We also point out how important features of the QCD axial anomaly 
are mirrored in this model. 

In the Hamiltonian approach, the above rf model has the following pairs of canonically 
conjugate fields: 

(r/,7r) and (B",P X ). (5) 
In addition, it has the first class constraints 

G(c) = J d 3 x (d lC l )(x) (P - At/) (x,t), (6) 

<?(„) = /«?,„«(*)«(*,,), (7) 

where c l ,w l are Schwartz functions. 

When there are Nf flavors, chiral symmetry breaking implies the existence of an Nf x Nf 
matrix of fields 

u = s e ir1 ', (8) 

where 

u(x) G U(N f ), s(x) G SU(Nf), e ir '\x) G 17(1). (9) 

The field s describes the (Nj — l)-dimensional flavor multiplet and e lTI ' describes the r/' field 
as in ([5]). 

The chiral group SU(Nf) L x SU(Nf) R acts on u according to the rule 

(9l,9r)>u(x) ^ g L u(x) g j w (10) 
(9l,9r) e SU(N f ) L x SU(N f ) R , (11) 
while the action of the axial vector group is 

u(x) —> u u(x), \oj\ = 1. (12) 
Now, the fields ^ are invariant under 

s(x) — >■ s(x) e N f , 

rf(x)^rf(x)-k—, (13) 
l\ f 

with k G {0, 1, 2, JV> — 1}. 
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That means in particular that the constraints and observables in field theory should be 



invariant under (13). The constraints (|6j |7j) already are so. 
Let D l be smooth functions on R 3 fulfilling the condition 



lim / dQ r 2 x.D^x) = 1. (14) 



Then 



Q(D) = j d 3 x (d^ix) (p - Ar/)(£,t) (15) 



commutes with the constraints. However it is not a constraint as D l are not all Schwartz 



functions. But it fails to be invariant under (13). However 



W(x) = s(x) e -^ Q{D) := s(x)V (16) 



is invariant under (13). So W(x) is an observable. 



We need to consider the following observables as well in what follows: 

a) J d?x B°(x,t) and its exponential form 

U{9) =e ^fd 3 *B°m. (17) 

b) 7T-XB . 

Let 

\-)e=o ■= |")o (18) 

denote any QCD #-state vector for cos# = 1 or 9 = mod 2it (Hereafter, whenever we 
write 6, we mean 9 mod 2tt.). Then the corresponding state vector for any cos# is 

U{9) |-)o. (19) 

We can write the fermion coupling to chiral fields either for quarks or for baryons. 
Let us first focus on quarks. The relevant couplings, chirally invariant, but with P and 
T violations are the quark mass terms: 

£q = Ml (ql u q R + q R u ] q L ) + (jl 2 (qi W q R + q R W ] q L ) . (20) 



The second term is new for 9 ^ 0. But for 9 = 0, W = u and (20) becomes the standard 
quark mass term with mass + /i 2 . 



Let 



be the axial charge of quarks: 



Qa = J d 3 x (q{q L - q R q R ) (21) 

e iaQA q L e~ iaQA = e~ ia q L , (22) 
e iaQA q R e~ iaQA = e +ia q R . (23) 



It is not a symmetry of (20). But just as in QCD we can have a "consistent" axial charge 



[SJ IS] by the addition to Qa the following term: 

Q w = -2 J d 3 x ir{x,t). (24) 

Then since 

e i*Q. e -iv' e -i*Q« = e 2ia e -ir,'^ (25) 

we have that 

Qa + Qw = constant of motion. (26) 
A. Perturbation Theory 

The QCD state vector |-)o is one which does not lead to P- or T-violation. Its vacuum 
vector |0)o is P- and T-invariant. If |-)o is a bipartite or multi-partite system with gluons, 
quarks and photons, neither the gluon nor the quark state vectors in |-) lead to P- or 
T-violation. 

In our model, the constituents in the multi-partite system are chiral fields and quarks. 
The dressing transformation U (6) does not contain quark fields. We will do perturbation 
theory in the electromagnetic coupling constant e. In that case, the free external quark state 
vectors are not twisted by U{6), only the chiral fields are. This point is important in what 
follows. 

In our calculation of EDM, we do not consider propagating chiral fields. It is enough thus 
to focus on the terms in the Hamiltonian involving the fermion fields. Their kinetic terms 
and coupling to electromagnetism have no chiral fields and are not affected by U{9). The 
terms in the Hamiltonian sensitive to U(8) are the quark mass terms: 

h° q = - J d 3 x [m (q L uq R + q R u ] q L ) + /i 2 (q L W q R + q R W ] q L )] (27) 



We now replace h° q by 

U(6)- 1 h° q U(6) = h e q = - J d 3 x [q L (m u + e~ ie fi 2 W) q R + h.c.] (28) 

in order that henceforth we can work with untwisted 8 = state vectors |-)o. 
In perturbation theory in e, to zeroth order, |-) splits as follows: 

|-)o = |chiral particles) o <8> | quarks and photons) o- (29) 

As we do not consider excitations of chiral particles, we can replace the first factor by 
|0)o, the vacuum of the chiral fields. In this vacuum, with the usual boundary conditions 

u(x,t) —7-1 as \x\ — > oo (30) 

we have 

o(OK£,i)|0) = <0|eW)|0)o = l. (31) 

In contrast to (48) of [1] which is appropriate for ^-states, we now have also 
o<0|P (aO|0) = 0. 

The expectation value of h q for the vector |0)o gives the following effective quark mass 
term H q with ^-effects included: 

H e q = o(0|/ij|0) = - J d 3 x [q L {fi 1 + e~ ie fi 2 )q R + h.c] 

= - J d 3 xq (/ix + e^V 2 ) q- (32) 



Set 



where 



Then 



0i + e*»V2 = M0)+H&/A0), ( 33 ) 

fi(9) = + fi 2 cos (34) 

H\Q) = / i 2 siii0. (35) 

Hx + e i75 V 2 = \M{9)\ e ^ m , (36) 
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with 



\M(8)\ 2 = f,(8) 2 + (,'(e) 2 , 
H{0) 



cosA(fl) 
sinA(fl) 



\M(9)[ 



(37) 
(38) 

(39) 



\M{9)[ 

The significance of \M{9)\ is that it is the quark mass renormalized by 9-angle. This can 



(40) 



be seen from the poles of the quark propagator for the mass term (33): 

-» 7 ■ k + M(0) 
F[ } fe2 + |M(0)|2-ze' 

where 

M{9) = \M{9)\ e ~^ m . (41) 

For the external spinors, it is thus natural to use the P- and T-invariant Dirac equation 
for mass \M{9)\: 

(ij- k+ \M(9)\) u(k) = 0. (42) 

We will do so hereafter. 



III. THE LOOP DIAGRAM 



The diagram we must calculate is a simple modification of the one which corrects the 
electron magnetic moment. It is shown in figure [TJ Following Weinberg |10] , section 11.3, 
page 485, we obtain the vertex function 



where 

rV,p) = / d 4 k (e 7 p (27r) 4 ) 



u{p') r M (p',p) u{p) 

i -17 • (p' -k) + M{9) 
(2tt) 4 ( P / - k) 2 + \M{9)\ 2 -ie 
i -«7 • (p - k) + M{9) 



(43) 



(44) 



(e7p(2vr) 4 ) 



(2tt) 4 k 2 - ie 



(2tt) 4 (p - A;) 2 + \M(9)\ 2 -ie_ 
This differs in an essential way from Weinberg's equation (11.3.1) only in that M{9) has 
a 75-term. We can thus simplify T^{p' ,p) exactly as in that book arriving at the analogue 
of his (11.3.4): 



(2vr) 4 7 J 



[k 2 + \M{9)\ 2 x 2 + q 2 y{x - y) - ie 



3 ' 



(45) 



where q = p' — p and 

f " =y [-i (7 ■ p'(i - y) - 7 • * - 7 • p(z - y)) + M{&)] 

Y [-i( 1 -p{l-x + y)- 1 -k- 1 -p'y) + M(e)\ lp . (46) 




Figure 1: Vertex correction diagram. 

The P- and T- violating term comes from the dependence of M{9) in on e - * 75 ^. Only 
the homogeneous, first order terms in M{6) retain these factors: they cancel in the second 
order one because of the intervening / y fl : 

M{6)YM{6) = \M{9)yf\M{9)\. (47) 



Retaining just first order M(#)-terms in (46) and dropping 7 ■ k terms ("symmetrical inte- 
gration"), we obtain T^y. 

f> -> f f 1} = -iM(6) 7 P 7" [7 • P (1 - x + y) - 7 • p' y] 1 P 

- iM(9) Y [7 • p' (1 - y) - 7 • p(a; - j/)] 7 M 7p (48) 

A. Averaging over x and y 



Following Weinberg [TU], suppose we change variables in the x,y integrations in (45) 
according to 

x 1 — y x, y x — y. (49) 



This does not affect the denominator in (45), nor does it affect the x-, y-integration in that 



equation. If we now add to (46) its expression after these changes and divide by 1/2, the 
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linear terms in T M simplify considerably letting us replace linear x- and y-terms as follows: 

x^x, y -(y + x-y) = |. (50) 



Of course, this rule is not valid for quadratic expressions in x and y, but by (48), we do not 
have to deal with such terms. 

Averaging over x and y simplifies to r£L: 



x 
2J 



1p 



— iM(6)^ p 7 ■ p — 7 • {p + p 
Since 7 P 7 M 7 A 7 P = 4^ , we can simplify this further: 



x 



(51) 



1 (2) 



We now use the identity 



U M{9) (p + p 



d^k (2tt)' 
(fc 2 + A) 3 ~ ~4A~ 



(52) 



(53) 



to obtain 



fi(p') ff 2) u(p) = 4e 2 tZ(p') M(9)(p+pY «(p) 

1 -z/2 



(54) 



/ y y |M(^)| 2 x 2 + g 2 ^-y)' 
The momentum dependence of this expression is governed by (p+p') M as required by current 
conservation. 



Note that we have not encountered ultraviolet divergence in arriving at (54). But the 



constant term in the numerator of the integrand of this integral, that is, the term 1, leads to 
infrared divergence. It is exactly the infrared divergence, which is homogeneous and linear 
in the electron mass, encountered in the anomalous magnetic moment calculation [11]. It 
is therefore exactly canceled by the contribution from soft photon emission processes |10j . 
The twist e 1 " 15 ^ in M{9) makes no difference. 



We can see this in another way. The result (54) can be interpreted as coming from 
changing incident and outgoing wave functions to 



u(p) 
u(p') 



ns x 2 e u{p), 
i75 ^ u(p'), 



(55) 
(56) 
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and working with |M(0)| and not with M(6). In this way the twist U(6) of if- (or "gluon- 
") vacuum has been transferred to w's. But this change does not affect the QED vertex: 
u(p')'y fJ- u(p) = u(p') , y tJ "u(p) . The cancellation of the infrared divergence in this transformed 
situation is explained in text books pU]. So it cancels here as well. 



Dropping the constant term in (54), we obtain the finite answer 



u(p') ff 2) u{p) -> u(p') ff 3l u(p) = -2e 2 u{p') \M{e)\e^ m u{p) (p + pY F{q 2 ), (57) 



(3) 



with 



F(q 2 



dx / dy 



x 



o J y \M(e)\W + q*y(x- y y 



Equation (57) has a P- and T-conserving part 



u(p') t^ E u(p) = -2e 2 |M(0)|cos(A(0)) u(p')u(p) (p+p'Y F(q 2 ) 



(58) 



(59) 



which has O (sin 2 6) corrections to the standard electromagnetic vertex. It has also a P- and 
T-violating EDM term 



u(p') rf 3)0 u(p) = -2e 2 |M(0)|sin(A(0)) u(p% l5 )u(p) (p + pY F(q 2 ). (60) 



Since 



F(0) 



we obtain the static electric dipole moment 

d = -2e 



|M(0)|* 



2 sin(A(0)) 



(61) 



(62) 



\M(9)\ ■ 

For the neutron, identifying |M(0)| with the neutron mass fi^, we obtain its EDM d^ as 



d N = -2e z 



sin(A(0)) 



— ona 



\i<i sin 8 

2 ' 
A* AT 



(63) 



where in the last equality we have used ( 39 ) and ( 35 ) . 



If we recall that experimentally |djv| < 6 • 10 26 cm, this gives the bound (for hn 
940 MeV « 7.5 ■ 10 12 cm" 1 , a = e 2 /47r « 1/137) 



sin(A(^))| = / ^^<2.5-10- 12 . 



(64) 
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IV. ON USE OF MIXED STATES 



We have recalled in Section II that if \-)e=o is any hadronic vector state for 9 = 0, then 
the corresponding vector state for finite 9 is 

\-)e = U(e)\-) e = . (65) 



The expression (63) for the EDM comes from the P- and T- violating part of the matrix 
element 

e(p\^l^\p)e- (66) 
Here we have suppressed labels like spin components and charge. 



The expression (66) is the expectation value of the current ^7^^ for the density matrix 



Pe = \p)e e(p\ ■ (67) 

e(p$YMp)e = Trpe^Y^ (68) 

(The spin sum on the RHS will introduce overall factors. We can ignore them.). 
The P- and T-transform of U(9) is U(—9). Hence that of pe is p_#. 
Our proposal in a previous work [lj was to use the mixed state 

* Pe + P-e , ra , 

P = g ( 69 ) 

as a possible mechanism to obtain zero EDM and solve the strong CP problem. 



We can now see that the proposal works. Since the matrix element (67) is proportional 
to sin#, it is clear that the EDM term in Tr p ip^^ip is zero. That is so no matter what sin# 
is. This means that with the use of the above mixed state, there is no way to constraint sin^ 
from EDM. There is no need for an axion either. 

Note that 

p^e ± PpeP ] = PpeP (70) 

That is because in p_#, we reverse only 9, whereas P will also change variables like momen- 
tum and helicity. For example, 

P\p)e e(p\P = \po, ~P)-e -e(po,p\, (71) 

whereas 

P-e = \p)-e -e(p\- (72) 
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Remarks 



Since \p)g is not normalizable, we should really work with 

/(M \k)e=\f)e, (73) 
where / is a function with compact support in k centered at p and normalized to 1: 

/£l/(MI' = l. C74) 

Let 

pe(f) = \f)e e(f\- (75) 

Then the mixed density matrix is 

m = nn+fM . (7 6 ) 

The EDM vanishes for p(f) as well: 

Terms odd under P and T in Tr p 9 (f) ipj^ip = 0. (77) 
Hence our previous conclusions are not affected. 

V. FINAL REMARKS 

There are good reasons to develop a quantum theory where certain symmetries do not 
become anomalous. Thus, consider the case dealt with in the present paper, namely parity 
and time- reversal anomaly in QCD when cos 8 ^ 1. There is no experimental basis for this 
anomaly. It has led to the suggestion that there is a particle called the axion. There is no 
evidence for the latter either. 

There are more of these examples. The following deserves special mention. The color 
group becomes anomalous in the presence of nonabelian magnetic monopoles [7J [121 113] - 
That too is the fate of mapping class groups ("large diffeos") in certain quantum gravity 
theories [16J. 

We may not always tolerate such anomalies. The use of mixed states is one way to recover 
them as symmetries. 

The entropy created by these mixed states is not always small. The entropy of p(f) is 
roughly log 2 and that is per an elementary quantum state. It can easily be of the order 

13 



of 5 or 10 per quantum state depending on the volume of the anomalous group. In QCD 
plasma, entropy per particle is about 10 |14j . The mixed state entropy can be comparable 
to 10. They can thus affect phenomena in elementary particle physics. 

Let <p n g) </?Jj be the resolution of identity corresponding to the eigenstates of the 
Hamiltonian for a particular domain T>. If a symmetry such as parity is anomalous, we 
should average this resolution over the anomalous group. For parity in our problem, the 
averaged resolution of identity is 

X - J2 {<P n ® & + P^n ® VIP' 1 ) ■ (78) 

n 

Call such an averaged resolution of identity in general as p. 

The use of anomaly-free mixed states implies that the Gibbs state is changed, the expec- 
tation values of correlators being 



Trp e~P H O x ---O n 
Tr p e~ 



(Oi-o N )= m/ ,,,, ^ w • (79) 



The different terms obtained by averaging in p can affect these expectation values as they 
come from different domains of the Hamiltonian H. 

This formula can also be adapted to thermofield theory [TS]. The properties of those 
correlators seem different from these for the standard Gibbs state. 

Nonabelian gauge symmetries G which define "charges" (as contrasted to those implied 
by Gauss law which hence vanish on vector states) have unusual properties [T7j. Observ- 
ables must commute with these symmetries. That means that only the elements in the 
center C(CG) of the group algebra CG are observable. By the preceding remark, if a gauge 
transformation g G G does not commute with C(CG), it is not observable. Restricting any 
pure state with nontrivial response to G to its center C(CG) creates a mixed state and en- 
tropy. We will develop this remark elsewhere. It already happens for the molecule C2-H4 
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